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Abstract. In radiative layers of rotating stars the luminosity carried by circulation currents through a surface of constant 
entropy (circulation luminosity) is shown to be positive. The corresponding decrease in the temperature gradient is important 
in the secondary of contact binaries. This result removes the deadlock in the theory of contact binaries. 

The resulting treatment of contact binaries is investigated, assuming thermal equilibrium. The sources of the circulation lu- 
minosity in the secondary can be written as the product of a circulation function (a normalized non-negative function of the 
fractional mass) and an amplitude. If the amplitude is adjusted to give a prescribed temperature difference ATc — Tci — Tc2, 
the choice of the circulation function is (in a broad range) unimportant. This invariance extends in a close approximation to 
all observable properties as well as to the internal structure. The temperature difference is bound to be positive. The fractional 
extent of radiative regions is larger in the secondary than in the primary. In the course of evolution the period increases and the 
mass ratio decreases. Comparing thermodynamic quantites on level surfaces, pressure and density are larger in the secondary 
than in the primary. The specific entropy is larger in the primary than in the secondary. The temperature difference is remarkably 
small and almost vanishing when averaged over the level surfaces occupied in both components. 

The only free parameter (apart from ATc) is the efficiency /e of the energy transfer from the primary to the secondary. 
Using standard values for the parameters, a survey of unevolved and evolved contact configurations is presented. Observational 
tests are passed. In stable systems the degree of contact is small. Stable systems in the period-colour diagram, unevolved 
and evolved, cover the strip (and only the strip) of observed systems in this diagram. Lower limits for period and effective 
temperature, compatible with the observed limits, are caused by the requirement of thermal stability. Stable systems with mass 
ratios very close to unity are possible, in accordance with recent observations. Since stability considerations are essential in 
these observational tests the results support the assumption of thermal equilibrium as well as the treatment of the stability 
problem. 

Models for individual observed systems with reliable data are well-determined (apart from some freedom in ATc) and can be 
used to calibrate the efficiency and to determine metallicity and age. All models obtained so far are stable. This again supports 
the assumption of thermal stability. The results show that evolutionary effects are important and that the efficiency is very small 
(/e = 10-^ . . . 10-^). 

Arguments are presented that the velocity field in the common envelope has a reversing layer, with motions from the secondary 
to the primary in the layers just above the critical surface and from the primary to the secondary in the surface layers. 

Key words, stars: binaries: close - stars: rotation 



1. Introduction 

This paper is the third one in a systematic attempt to overcome 
the deadlock in the theoretial treatment of contact binaries. In 
the first paper (Kahler 120021 hereafter Kl) we discussed the 
uncertainties in the structure equations. Serious uncertainties 
were found to concern only the energy sources and sinks asso- 
ciated with internal mass motions. In the second paper (Kahler 
12002 ' hereafter K2) we made the usual assumption (hereafter 
assumption A) that sources and sinks occur only in the com- 
mon envelope. Investigating a typical system we found a va- 
riety of solutions which are reasonable from an observational 
point of view, but all of them suffer from theoretical deficien- 
cies. In particular, all solutions with good light curves (i.e. a 
small temperature difference between the components) are of 
the type proposed by Moss & Whelan (il973< . i.e. the energy 



sources in the secondary's envelope occur in the extreme su- 
peradiabatic layers. As shown by Hazlehurst J1974> this is im- 
possible since the heat capacity of these layers is too small. 

The results in K2 show that assumption A is too restric- 
tive. When thermal equilibrium is assumed this conclusion is 
cogent since it does not depend on details, e.g. the choice of a 
transport equation for the energy transfer from the primary to 
the secondary. In the case of thermal cycles there is more free- 
dom since they depend on the transport equation as well as on 
the equation governing the rate of mass exchange. Nevertheless 
the experience of decades shows that the conclusion is justified 
also in this case. Associated to a cycle is a configuration in un- 
stable thermal equilibrium. If this configuration is in shallow 
contact (and if extreme superadiabatic transfer is excluded) the 
temperature difference is large for a considerable part of the 
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cycle. If the configuration is in good contact the temperature 
difference is always large. In both cases the light curve is bad. 

Assumption A turns out to be a deadlock in the treatment 
of contact binaries. Energy sources or sinks in the interior of at 
least one component and corresponding changes in the temper- 
ature gradient are necessary. According to Lucy (J968 l there 
is no freedom in the structure of radiative regions since the lu- 
minosity carried by circulation currents through a level surface 
vanishes (hereafter Lucy's theorem). The only remaining free- 
dom concerns the extent of turbulent regions. A decrease of a 
convective envelope is in conflict with Lucy's theorem. An in- 
crease of the turbulent envelope in the secondary is in conflict 
with Roche geometry. An increase in the primary is compati- 
ble with Roche geometry, but the resulting model (Kahler l 19891 
Kahler & Fehlberg 1991 1 is in conflict with the second law of 
fliermodynamics (Hazlehurst ll993> since (in a large region) in- 
ternal mass motions transfer energy from cooler to hotter lay- 
ers. 

Accordingly, the deadlock boils down to Lucy's theorem. 
This theorem has been derived — also by Mestel J1965> and 
Roxburgh et al. (1965 1 — assuming that all thermodynamic 
quantities are constant on equipotential surfaces, i.e. assum- 
ing pseudo-barotropic layers in strict hydrostatic equilibrium. 
Although the theorem is usually taken for granted (e.g. Maeder 
& Zahn .1998) . it is valid only in the zero-order approximation 
and not in higher approximations (Mestel 1965 ). In contact bi- 
naries departures from strict hydrostatic equilibrium (and thus 
baroclinic effects) are necessary for the energy transfer from 
the primary to the secondary, and baroclinic effects are ex- 
pected to be important also in the interior of the components 
since the variation of the gravity on a level surface is large. For 
these reasons it is likely that Lucy's theorem does not apply. 
As observed also by Tassoul & Tassoul ('I995'l in a discussion 
of rotating stars there is no reason to claim that the net flux 
of energy carried by the circulation currents through each level 
surface vanishes. 

Removing the deadlock and assuming thermal equilib- 
rium we are able to predict the modification of the original 
Lucy {1968 1 model which solves the contact binary problem. 
Roche geometry requires that (in comparison with the Lucy 
model) either the (averaged) temperature gradient in the pri- 
mary be increased or that the gradient in the secondary be de- 
creased. An enhanced gradient is in conflict with the second 
law. Accordingly we need a reduced gradient in the secondary 
and thus a reduced radiative luminosity and/or a reduced extent 
of convective layers. The reduction is possible if some part of 
the luminosity in the secondary's interior is carried by circula- 
tion currents to the common envelope. 

In this paper we investigate the resulting treatment of con- 
tact binaries, assuming that thermal equilibrium is enabled 
by circulation currents in the secondary's interior We discuss 
the effects of circulation, present a survey of unevolved and 
evolved theoretical configurations, and perform observational 
tests. Having obtained models for individual observed systems 
we shall finally return to the hydrodynamic problem. 
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Fig. 1. The degree of contact F in dependence on the ampli- 
tude / of the circulation in System 1, for different values of the 
parameter a (see text). 
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Fig. 2. The temperature difference ATc in dependence on the 
amplitude / in System 1, for different values of the parameter 
a. 
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Fig. 3. ATe as a function of / for a = 0.95 in System 1. 
Branches 1 and 3 (solid lines) of the linear series are stable, 
while branch 2 (dashed line) is unstable (see text). 
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2. The effects of circulation currents 

We shall adhere to the notation of Kl and K2. Roche geom- 
etry is vahd in a close approximation and will be assumed 
throughout this paper, apart from the discussion of the veloc- 
ity field of the internal mass motions in Sect. 15.11 An equa- 
tion for the transfer of mass between the components will not 
be needed. The energy transfer from the primary to the sec- 
ondary will be treated as proposed in Kl, assuming ai =012- 
Recall that a denotes the extent of the energy sinks/sources 
in the primary's/secondary's envelope, expressed as fractional 
mass of the layers above the critical surface. Since extreme 
superadiabatic transfer is not possible, the case of a ^ 1 is 
excluded. Since the transfer occurs in the common envelope, 
values of a much larger than unity are also excluded. The pa- 
rameter /e, measuring the efficiency of the energy transfer, is 
certainly smaller than unity and possibly much smaller The 
mixing length in units of the pressure scale height will taken 
to be a+ = 1.5. We shall assume that the extent of turbulent 
regions is determined by the Schwarzschild criterion {it = 0). 

With these assumptions we have so far only two free param- 
eters a, /e. An additional parameter will be required when in- 
troducing circulation currents in the secondary. A modification 
of the structure equations derived in Kl will also be required. 

2. 1. Energy balance and temperature gradient 

In the common envelope the luminosity A is transferred from 
the primary to the secondary (hereafter process 1). Let us as- 
sume that in addition in the secondary a fraction / of A is trans- 
ferred by circulation currents from the interior to the envelope 
(process 2). Here we collect the equations describing the super- 
position of these processes. 

The local balance of energy in the component i is 



dli/ drrii 



(1) 



where cr, represents the sources or sinks in the envelope caused 
by process 1, i.e.by the energy transfer form the primary to the 
secondary in a region described by the parameter a. (Although 
thermal equilibrium is assumed, the time-dependent form of 
the energy balance is needed when discussing thermal stabil- 
ity.) Note that <J2 does not allow for process 2. 

In the presence of circulation in a radiative region of the 
secondary the luminosity is the sum 



rad,2 



+ ^circ,2 



(2) 



of the radiative luminosity ^iad,2 and the luminosity carried 
by circulation currents lciic.2 (hereafter circulation luminosity). 
The temperature gradient in a radiative layer 



'rad.2 



167racG ni2T2 



(3) 



is coupled with /iad.2 and not with ^2- (Thus, in the presence of 
circulation Eq. (5) in Kl is not valid.) The borders of convec- 
tive regions are therefore also coupled with ^rad,2- Cumulation 
occurs as well in convective regions (e.g. Kippenhahn I 1 963> . 

To investigate the effects of process 2 we shall assume a run 
of Zcirc,2 throughout the secondary as a non-negative function 



of the mass variable which vanishes at the boundaries. (Since 
process 2 amounts to a redistribution of energy, the sources 
of the circulation luminosity, i.e. the regions with a positive 
derivative d^circ,2/d™2 > 0, are compensated by sinks with a 
negative derivative.) Equations Q and (|3ji define the radiative 
gradient Viad,2- If this gradient is not larger than the adiabatic 
gradient the layer is radiative. Otherwise the layer is convective 
and the temperature gradient can be obtained from the mixing- 
length theory. Accordingly, the temperature gradient is deter- 
mined as usual, replacing simply the luminosity I2 by the ra- 
diative luminosity lrad,2- We have tacitly assumed that /rad.2 is 
positive apart from the centre. Otherwise the choice of /circ.2 is 
unphysical. 

Writing the circulation luminosity as a function of the frac- 
tional mass X2 = m2/M2, we shall assume that the derivative 
is a superposition 



d/circ,2/dX2 = /A 0(2:2) + 



(4) 



of sources and sinks. In the first term, representing the sources, 
c{x2) is a normalized non-negative function which will be 
called circulation function. The amplitude / is expected to be 
small compared to unity in realistic models with typical mass 
ratios. The dots represent corresponding sinks in the envelope 
with a distribution described again by the parameter a. 

2.2. The choice of the circulation function 

Since so far nothing is known about realistic circulation func- 
tions, two simple functions depending on a parameter a will be 
used. Writing for simplicity x instead of X2, the first function 



Cl 



{1 + cos[tt{x — a) /d]} / {2d) if x>a — d 
otherwise 



(5) 



with d ~ X{ — a (where X{ represents the fitting mass) has its 
maximum at a; = a. The second function 





C2 = •\ V {I - COS [7r(x - a)/{b - a)]} 
with 



if a; < a 

if a < X <b (6) 
ifx>b 



b={a + Xf)/2, V =l/(2xf-a-b) 



(7) 



has a non-negative derivative, vanishes at the centre, and in- 
creases between x — a and x — b. 

As in K2 we investigate an unevolved system with the phys- 
ical parameters 



M ^ 1.6M0, J52 = 0.4, X ^ 0.7, Z = 0.02, 



(8) 



where J52 denotes the angular momentum in 10 cgs-units 
(hereafter system 1). Using the parameters a — 0.5, /e = 0.1 
and the circulation functions ci(a, x) with different values for 
a we calculated a continuous sequence of configurations in de- 
pendence on the amplitude / of the circulation, i.e. a linear 
series with the parameter /. Results for the temperature differ- 
ence ATe = Tci—Tc2 and the degree of contact F are shown in 
Figs-IDandlJl In the absence of circulation ( / = 0) the tempera- 
ture difference is very large and the contact is shallow. We need 



H. Kahler: The structure of contact binaries 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 





■ ', 


\ 


; ! 


\ \ 
\ 

\ 


\ 0.1477 




\ 0.4 










0.8', 











0.15 
F 



Fig. 4. The temperature difference ATc as a function of F in Fig. 6. Lines of constant angular momentum in the (ATc, F) 
System 1, for the same values of the parameter a as in Figs.^ diagram for systems with M = I.GMq. The numbers denote 



andH 



the mass ratio for AT^ = 200 K (see text). 
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Fig. 5. ATc as a function of F in System 1, for a — 1 and 
different values of /e (solid lines), for a = 0.5, /e = 0.1 
(dotted line), and for a = 5, Je — 0.1 (dashed hne). 

configurations in shallow contact with a small temperature dif- 
ference. For any choice of a they can be obtained with a small 
amplitude (/ ~ 0.04 . . . 0.08). This is as expected from the ar- 
guments collected in Sect. \l\ Unfortunately the curves depend 
sensitively on a, i.e. on the choice of the circulation function 
which so far is unknown. 

In particular, note that the linear series has turning points 
if a is sufficiently large (a > 0.9), i.e. if sources of the 
circulation luminosity occur only in outer layers. For a = 
0.95 this is illustrated in Fig. |3] The curve consists of three 
branches. Configurations on branch 2 (dashed line) are unsta- 
ble. Configurations which are reasonable from an observational 
viewpoint (in the neighbourhood of the turning point between 
branches 2 and 3) are either unstable or close to instability. 
Unstable configurations are unrealistic. Solutions which are 
stable but close to instability are extremely fragile from a com- 
putational viewpoint and thus probably also unrealistic. This 
suggests that the sources of the circulation luminosity occur not 
only in outer layers but also deep in the secondary's interior. 



The uncertainty in the choice of the circulation function 
seems to be a serious drawback. However, if / is eliminated and 
ATc considered as a function of F, all curves coincide apart 
from very large values for F which are not of interest. This 
surprising result is shown in Fig.|3 In other words, the relation 
between the two observable quantities ATg and F does not de- 
pend on the choice of the circulation function. This invariance 
(which has been checked in numerous cases and which extends 
much further as we shall see) enables a treatment of contact 
binaries although so far we are not able to determine details of 
the circulation. 

Figure [S] shows the relation between ATc and F for sys- 
tem 1 in the region of interest (a small temperature difference 
and shallow contact) in more detail and for different values of 
the parameters a and /e. The dependence on a is weak and can 
be neglected in an approximate treatment, but the dependence 
on /e is important. The curve for /e = 1 represents an unreal- 
istic limiting case. Realistic values for /e are certainly smaller 
than unity and probably much smaller. The curves show there- 
fore that the temperature difference in system 1 is positive. 

2.3. The temperature difference 

Having seen that the temperature difference in system 1 is 
bound to be positive, here we ask whether this result can be 
generalized. As an example we consider unevolved systems of 
the same mass and composition as before, with a ~ 1, /e = 
0.1. Figure|6lshows the invariant relation between ATc and F 
in systems with different angular momentum. The curves are 
labeled by the mass ratio for ATc = 200 K. The curve begin- 
ning with ATc = 200 Kandq ^ 0.8 ends up with ATc ^ 
and q ^ 1, i.e. with a system with equal components. 

It is manifest that the temperature difference is positive 
apart from the limiting case q ^ 1. This result has been 
checked in many other cases. It is valid for all stable systems, 
evolved or unevolved. 

Concering a physically plausible interpretation, circulation 
currents tend to decrease the departures from hydrostatic equi- 
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Fig. 7. Luminosities (in solar units) in the secondary's interior 
in System 2 in dependence on the fractional mass. The stan- 
dard circulation function is used. The total luminosity (solid 
line) is the sum of the radiative luminosity (dotted line) and the 
circulation luminosity (dashed line). 

librium in the whole system. Therefore they tend to reduce the 
temperature difference, which means that 

dATe/d/ < 0, (9) 

as it is the case in the relevant parts of the curves in Fig.Qunless 
the branch is unstable, i.e. unless the parameter a is too large. 
Since the circulation cuiTents do not quite succeed in removing 
the temperature difference, this difference remains positive. 

In the example just considered the system with equal com- 
ponents is unstable, but systems with mass ratios up to g = 
0.9918 are stable. For a very long time the absence of systems 
with nearly equal components seemed to be an established ob- 
servational fact, and an instability of systems with q ~ 1 was 
discussed as a possible theoretical explanation. The present ex- 
ample shows that the existence of systems with mass ratios very 
close to unity cannot be excluded from a theoretical viewpoint. 
This is in accordance with the observation of the system V 753 
Mon with q ~ 0.97 (Rucinski et al. .2000^ 

2.4. The treatment of contact binaries 

The discussion in the preceding subsections has shown that in 
systems with typical mass ratios a small circulation luminosity 
(few percent of A) in the secondary's interior has a large influ- 
ence on the structure. A similar circulation luminosity in the 
primary's interior has almost no effect since the radiative lumi- 
nosity is much larger. For this reason circulation in the primary 
will be neglected, keeping in mind that in systems with large 
mass ratios the situation may be more complex. 

To allow for circulation in the secondary we have to spec- 
ify (i) a circulation function and (ii) the amplitude / or another 
measure for the strength of the circulation. Figure|2]shows that 
the amplitude which is necessary for a given temperature dif- 
ference depends on the choice of the circulation function. To 
obtain a measure which is invariant the parameter / will be 
replaced by ATe. In other words, the amphtude / will be ad- 
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Fig. 8. The circulation luminosity in the interior of the sec- 
ondary in System 2, in dependence on the fractional mass, for 
different circulation functions (see text). 
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Fig. 9. The entropy distribution in System 2 in dependence on 
the fractional radius, for different circulation functions. The 
upper curves, representing the primary, coincide. The lower 
curves represent the secondary. Flat parts indicate convective 
regions. 

justed to give the precsribed value for ATc. The procedure is 
unique since in the case of a reasonable circulation function 
(a < 0.85 in the example considered) and a positive tempera- 
ture difference the connection between amplitude and temper- 
ature difference is unique. 

The resulting treatment of contact binaries is essentially 
based on the assumption of thermal stability. For given physical 
parameters (and mixing length) the structure of a configuration 
depends on the choice of the free parameters a, /e, ATo and 
on the choice of the circulation function. How large are the re- 
sulting uncertainties? We begin using the values 

a = 1, /e 0.1, ATc = 200K. (10) 

Keeping the parameters fixed, let us ask whether the observ- 
able properties depend on the choice of the circulation function. 
Results for system 1 are listed in Tabled The first column con- 
tains the circulation function, the second one the parameter a. 
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Table 1. Observable properties of system 1 for a = 1, /e = 0.1, ATc = 200 K and different circulation functions (see text). 





a 


P 


q 




L 


Li 


L2 


Ri 


R2 


F 


Cl 


0.5 


0.2895 


0.4195 


5465 


1.125 


0.812 


0.313 


0.983 


0.657 


0.0341 


Cl 


0.6 


0.2890 


0.4196 


5462 


1.120 


0.808 


0.312 


0.982 


0.656 


0.0342 


Cl 


0.7 


0.2886 


0.4195 


5460 


1.116 


0.805 


0.311 


0.981 


0.656 


0.0342 


Cl 


0.8 


0.2883 


0.4196 


5459 


1.114 


0.804 


0.310 


0.980 


0.655 


0.0343 


C2 





0.2889 


0.4194 


5462 


1.119 


0.808 


0.311 


0.982 


0.656 


0.0342 


C2 


0.1 


0.2888 


0.4193 


5461 


1.118 


0.807 


0.311 


0.981 


0.656 


0.0343 


C2 


0.2 


0.2886 


0.4193 


5461 


1.117 


0.806 


0.311 


0.981 


0.656 


0.0343 


C2 


0.3 


0.2885 


0.4192 


5460 


1.116 


0.805 


0.310 


0.981 


0.656 


0.0343 



Table 2. Observable properties of system 2 for a = 1, /e = 0.1, ATo = 200 K and different circulation functions. 





a 


P 


q 




L 


ii 


L2 


Ri 


R2 


F 


Cl 


0.5 


0.4182 


0.3606 


6368 


4.024 


2.994 


1.030 


1.397 


0.873 


0.0599 


Cl 


0.6 


0.4179 


0.3606 


6367 


4.018 


2.989 


1.029 


1.396 


0.873 


0.0600 


Cl 


0.7 


0.4177 


0.3605 


6366 


4.014 


2.987 


1.028 


1.396 


0.872 


0.0600 


Cl 


0.8 


0.4176 


0.3603 


6366 


4.012 


2.985 


1.027 


1.396 


0.872 


0.0601 


C2 





0.4177 


0.3606 


6367 


4.017 


2.989 


1.029 


1.396 


0.872 


0.0599 


C2 


0.2 


0.4177 


0.3605 


6367 


4.015 


2.987 


1.028 


1.396 


0.872 


0.0599 


C2 


0.4 


0.4175 


0.3604 


6367 


4.013 


2.985 


1.027 


1.395 


0.872 


0.0600 


C2 


0.6 


0.4175 


0.3603 


6366 


4.012 


2.985 


1.026 


1.395 


0.872 


0.0601 



Table 3. Observable properties of system 1 in dependence on the free parameters. 



a 


/e 




P 


q 




L 


Li 


L2 


Ri 


R2 


F 


0.5 


0.1 


200 


0.290 


0.419 


5461 


1.122 


0.810 


0.312 


0.984 


0.657 


0.033 


1 


0.1 


200 


0.289 


0.419 


5462 


1.119 


0.808 


0.311 


0.982 


0.656 


0.034 


2 


0.1 


200 


0.288 


0.420 


5463 


1.118 


0.806 


0.311 


0.980 


0.656 


0.035 




0.0001 


200 


0.280 


0.428 


5435 


1.086 


0.776 


0.309 


0.971 


0.660 


0.141 




0.001 


200 


0.284 


0.419 


5450 


1.101 


0.791 


0.310 


0.975 


0.658 


0.088 




0.01 


200 


0.287 


0.419 


5458 


1.111 


0.800 


0.311 


0.978 


0.656 


0.055 




0.1 


200 


0.289 


0.420 


5462 


1.119 


0.808 


0.311 


0.982 


0.656 


0.034 




0.1 


100 


0.285 


0.423 


5469 


1.106 


0.779 


0.327 


0.972 


0.653 


0.040 




0.1 


200 


0.289 


0.419 


5462 


1.119 


0.808 


0.311 


0.982 


0.656 


0.034 




0.1 


300 


0.292 


0.416 


5455 


1.131 


0.835 


0.296 


0.990 


0.659 


0.031 



Listed are the period P (in days), the mass ratio q, the mean 
effective temperature To, the total luminosity L (in solar units), 
the luminosities and radii of the components (in solar units), 
and the degree of contact F. Table |2l contains the correspond- 
ing results for the system with 

M = 2Mq, J52 = 0.6, X = 0.7, Z = 0.02, (11) 

(hereafter system 2). Larger values for a than those included 
in the tables lead to configurations which are either close to 
instability or unstable. This shows again that the sources of the 
circulation extend deep into the secondary's interior. 

Tables^andl^show that the results for different circulation 
functions are in close agreement. The function ci with a = 0.5 



represents a limiting case (a has the smallest possible value). If 
the results for this function (first line in the tables) is omitted, 
the agreement becomes excellent. We conclude that the choice 
of the circulation function (in a broad range) is unimportant as 
far as observable quantities are concerned. In other words, the 
invariance encountered in Sect. l2.2l extends in a close approx- 
imation to all observable properties. This invariance suggests 
the use of a standard circulation function. We decided to use 
the function C2 with a — 0. 

Consider next the effects of uncertainties in the free param- 
eters. Results for system 1, obtained using the standard circula- 
tion function, are listed in Tabled which is expected to cover 
the range of possible variations of the parameters. Changes in 
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Fig. 10. The log P — log T diagram in System 2, for different 
circulation functions. P is the pressure in cgs-units. The upper 
curves, representing the primary, coincide. The lower curves 
represent the secondary. 
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Fig. 11. logP as a function of the potential in System 2, for 
different circulation functions (see text). The lower curves, rep- 
resenting the primary, coincide. The upper curves represent the 
secondary. 



a have little influence, in accordance with a result from Fig.|5] 
An decrease in the efficiency /e leads mainly to an increase 
in F (as it might have been expected) and to a decrease in the 
primary's luminosity. A change in the temperature difference 
affects aU observable properties, particularly the light ratio. 

Since these effects are only moderate, an approximate treat- 
ment of contact binaries is possible using standard values of the 
free parameters. When getting a first survey of contact configu- 
rations in Sect.|3lwe shall use the values given in Eq. (I10> . In a 
next step observational tests can be used to improve these stan- 
dard values. In particular, Table|3lshows that reliable observed 
values for F can be used to calibrate the efficiency. This will 
be done in Sect.E] 
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Fig. 12. log T as a function of the potential in System 2, for dif- 
ferent circulation functions. The solid curves, representing the 
primary, coincide. The dotted curves represent the secondary. 
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Fig. 13. The specific entropy s as a function of the potential in 
System 2, for different circulation functions. The sohd curves, 
representing the primary, coincide. The dotted curves represent 
the secondary. 



2.5. The internal structure of the components 

We shall first consider system 2 as a typical example for the 
internal structure of the components and for the influence of 
the choice of the circulation function. Figure Q shows the lu- 
minosity in the secondary's interior as the sum of the radiative 
luminosity and the circulation luminosity. (The increase of the 
luminosity in the envelope escapes notice in this diagram.) The 
standard circulation function is used. Figure |8] shows the circu- 
lation luminosity in the secondary's interior for all circulation 
functions listed in Table |2] It is manifest that two classes of 
functions are used. The curve for the function ci with a = 0.5 
(solid line) is untypical since the amplitude / is rather large. 
As mentioned already this function represents a limiting case 
which is unlikely to be realistic. In the further discussion this 
function will be omitted. In the following diagrams for sys- 
tem 2 the remaining seven circulation functions are used. 
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Fig. 14. The mass ratio q in dependence on the angular mo- 
mentum J52 for unevolved systems with Z = 0.02 (see text). 
The numbers denote the mass in solar units. Solid/dahed lines 
represent systems of positive/negative charge. 

Figure |9] shows the distribution of entropy (per nucleus di- 
vided by Boltzmann's constant). Flat parts of the curves indi- 
cate convective regions. The curves for the primary coincide. 
The curves for the secondary coincide only in the envelope. 
The uncertainty in deeper layers is not large. Note that the frac- 
tional extent in radius of the convective envelope is larger in 
the primary than in the secondary. 

The logP — logT diagram (Fig. \10i illustrates the close 
correlation of the two components in the outer layers. Again 
the curves for the primary coincide. The uncertainties in the 
secondary's curves are remarkably small. 

Recall that in both components an effective central poten- 
tial is well-defined, c.f. Kl . Near the surface this potential coin- 
cides with the spherically averaged Roche potential. Recall fur- 
thermore that in the idealized case of strict hydrostatic equilib- 
rium all thermodynamic quantities are constant on connected 
equipotential surfaces. Departures from strict equilibrium are 
necessary for the energy transfer In particular, it is possible 
that circulation currents in the common envelope are driven by 
pressure differences on equipotential surfaces. Concerning the 
adjustment of circulation currents in the interior energy consid- 
erations might be important. 

For these reasons it is of interest to investigate thermody- 
namic quantities in dependence of the potential. We shall use 
the normalized potential which is positive, somewhat smaller 
than 2 at the surface, and larger in the interior. FigurelTTIshows 
logP as a function of the potential. The pressure difference 
AP = Pi — P2 is negative throughout the system, including 
the outer layers. We found that this result is valid for any sys- 
tem, unevolved or evolved, for any mass and any mass ratio. 
Implications for the mass motions in the envelope will be dis- 
cussed in Sect. 15. II 

Figure^lshows log T in dependence on the potential. The 
temperature difference AT ^ Ti — T2 is positive in the outer- 
most layers, changes sign in deeper layers, and is always small. 
When averaged over all layers the difference is very small. This 
surprising result was again found to be valid for any system. 



Fig. 15. The period-mass ratio diagram for unevolved systems 
with Z = 0.02. The symbols are as in Fig. 1141 




Fig. 16. The mean effective temperature in dependence on the 
mass ratio q for unevolved systems with Z = 0.02. 



Combining the results for AP and AT and making use of 
the equation of state we conclude that the difference in den- 
sity Ap — pi — p2 is negative throughout the system. Consider 
finally the entropy difference As = si — S2 (Fig. \131 . This 
difference is positive everywhere and large apart from the en- 
velope. 

From these results concerning differences of thermody- 
namic quantities on equipotential surfaces we see that most 
differences are large and have a characteristic sign. An excep- 
tion is the difference in temperature which is small and almost 
vanishing when averaged over the level surfaces occupied in 
both components. The temperature is closely connected with 
the internal energy. A connection between the total energy of 
the system and the thermal adjustment suggests itself. We shall 
return to this point in a forthcoming paper. 

This concludes the discussion of system 2. Results for sys- 
tem 1 are similar, with slightly larger uncertainties. Both sys- 
tems have intermediate masses and mass ratios, and both sys- 
tems are stable. When passing over to systems with smaller 
masses and/or smaller mass ratios, the uncertainties may in- 
crease and the stability may be lost. An example is the system 
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Fig. 17. logL (where L is the total luminosity in solar units) Fig. 19. Unevolved systems with Z = 0.02 (heavy Unes) and 
in dependence on the mass ratio q for unevolved systems with Z ~ 0.01 (weak lines) in the period-mass ratio diagram. 



Z = 0.02. 



Solid/dashed lines indicate configurations of positive/negative 
charge. The numbers denote the mass in solar units. 




0-2 0.4 0.6 0.8 1 

q 

Fig. 18. The degree of contact F in dependence on the mass 
ratio q for unevolved systems with Z = 0.02. 

with M = IMq and q = 0.32. Since the secondary's mass is 
very small there is a strong tendency towards convection and 
thus little freedom for the effects of circulation in radiative re- 
gions. Entropy distribution and number and extent of convec- 
tive regions depend strongly on the choice of the circulation 
function. Accordingly, the uncertainties in the internal structure 
are rather large in this case. These uncertainties do not extend 
to observable quantities. They are well defined. Uncertainties 
in the interior are not of interest since all configurations are 
unstable. 

Summarizing, in stable systems the invariance found in the 
preceding subsections extends also to the internal structure of 
the components. 

2.6. Remarks on the thermal stability problem 

Recall first how the thermal stability of a configuration can be 
tested. A full test requires either a stability analysis (i.e. the 
inspection of eigenvalues) or evolutionary calculations, start- 
ing with small perturbations. A restricted test requires only the 



charge c of a configuration, which is the sign of the Henyey 
determinant apart from a factor depending on the choice of the 
Henyey matrix. The charge c = 1 is necessary (but not suffi- 
cient) for stability, and the charge c = — 1 indicates instability. 

Consider the second test, i.e. the evolution following a 
small perturbation. The evolution can be calculated if the bal- 
ance of energy and the temperature gradient are known at each 
point of time. In other words, we need the sources/sinks ai and 
the run of the circulation luminosity at each point of time. 

The sources/sinks a-i are approximately known if the trans- 
port equation for A is reasonable. Assuming that the transport 
equation here used [Eq. (49) in Kl] is reasonable, this means 
that the choice of /e is reasonable. The circulation luminosity 
is known if the circulation function (represented by the param- 
eter a) and the amplitude / are known. Since there is a large 
freedom in the choice of the circulation function a can be kept 
fixed, but changes in / may be important. 

Throughout this paper we assume that / is kept fixed during 
a perturbation. The resulting stability problem is well-defined 
but simplified. To make this point clear note that the present 
definition of the amplitude / is arbitrary and that other def- 
initions are possible. For example, the expression /A in the 
RHS of Eq. (0} could be replaced by fL or by /Lq or by 
some other expression. All these possibilities are equivalent as 
far as the structure of a configuration in thermal equilibrium 
is concerned since the amplitude is adjusted to the prescribed 
temperature difference. They are however not equivalent in the 
stability problem with a fixed amplitude. The resulting uncer- 
tainty in the stability problem remains to be investigated, but 
the results obtained so far suggest that the stability problem 
here used is reasonable. 

All configurations of positive charge which have been 
tested so far turned out to be stable. For simplicity a config- 
uration is therefore called stable if c = 1. If the stability has 
been tested by evolutionary calculations, this is mentioned. 
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Fig. 20. Evolutionary effects in the period-mass ratio diagram 
for systems with M = IMq and Z = 0.02 (see text). Shown 
are curves of constant hydrogen content in the primary's centre 
(heavy lines) and curves of constant angular momentum (dotted 
lines). The numbers denote Xd and J52, respectively. Solid 
and dashed lines are defined as in 
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Fig. 21. Same as Fig.|20l but for systems with M = I.SMq. 

3. A survey of contact configurations 

In this section we present a survey of contact configurations, 
adopting the standard circulation function and the standard val- 
ues ( I10> for the free parameters. Under these circumstances an 
unevolved configuration depends on three parameters (mass, 
angular momentum, metallicity). An evolved configuration de- 
pends in addition on the chemical profile in both components. 
In the simplest treatment of systems with intermediate or small 
mass ratios a simple chemical profile is assumed and evolu- 
tionary effects in the secondary are neglected. The evolution- 
ary status is then described by the primary's central hydrogen 
content, and we are left with four parameters. 

Since a complete survey is impossible we decided to dis- 
cuss all unevolved configurations (in a certain mass range) 
with a standard composition, and to illustrate evolutionary ef- 
fects and the effects of changes in metallicity in few examples. 
Further examples are given by the observed systems discussed 



in Sect. |4] Examples for the effects of changes in the free pa- 
rameters have already been discussed. 

3.1. Unevolved configurations witti Z = 0.02 

We investigate systems with X = 0.07, Z = 0.02 in the mass 
range 1 < M/Mq < 2.5. The lower mass limit is caused 
by the lack of opacities for low temperatures and large den- 
sities. [The opacities (OPAL 95 combined with opacities of 
Dave Alexander) have kindly been provided by A. Weiss.] For 
M < 1.5Mq the opacities have been slightly extrapolated. 
For this reason and since complications in the equation of state 
(e.g. Coulomb interactions) are ignored, the treatment of sys- 
tems with small masses is approximate only. 

Figure [141 shows the mass ratio g as a function of the an- 
gular momentum J52 for different masses. The upper border 
of the region covered by configurations (near q = 0.9) is de- 
termined by the condition that the amplitude of the circulation 
is non-negative. In other words, the upper border is the curve 
f — which is invariant. For smaller values of the temperature 
difference the upper border occurs at larger mass ratios. 

The lower border is the line g = 0.1 since the Chebyshev 
approximation for the neck used in the numerical code (c.f. Kl) 
is valid only for q > 0.1. An exception concerns small masses 
(M < 1.37\/q). For these masses the lower border, occuring at 
somewhat larger mass ratios, is caused by the prescribed tem- 
perature difference. Below the border the difference is bound 
to be larger than 200 K. 

The relation between angular momentum and mass ratio for 
given mass is unique. Since the derivative dq/dJ52 is positive, 
the system with the largest mass ratio has the largest angular 
momentum. Systems of sufficiently large angular momentum 
(or mass ratio) are stable, while the remaining systems are un- 
stable. 

When passing over from intermediate to small masses the 
border between stable and unstable systems is shifting towards 
large mass ratios. For a mass somewhat smaller than the so- 
lar mass stability becomes impossible. This holds as well for 
smaller values of AT^. Accordingly, there is a lower limit for 
the mass of unevolved contact binaries, and thus also a lower 
limit for the effective temperature. The limit is caused not by 
the lack of configurations but by their instability. 

Figure[2]shows the period-mass ratio diagram. Solid lines 
represent again systems of positive charge. They do not cross, 
which implies that the structure of an unevolved system in ther- 
mal equilibrium with Z = 0.02 is (in a reasonable approxi- 
mation) fully determined by two observable properties. In the 
mass range under consideration there is a minimum period of 
about 0.18 days which is compatible with the observed mini- 
mum of 0.221 days for CC Com. The maximum period (about 
0.48 days) however is far below the observed maximum. 

Figures 11611181 show characteristic properties in depen- 
dence on mass and mass ratio. The results for the degree of 
contact F offer an observational test. Observed values for F 
are usually small compared to unity. This is in accordance with 
the results in Fig. ^] since for all stable systems F is small 
compared to unity. This lends support for the treatment of the 
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Fig. 22. The period-colour relation for unevolved systems with 
Z = 0.02. Shown are lines of constant mass (strong lines, 
solid/dashed for positive/negative charge) and of constant mass 
ration (weak lines). Observed systems occur in the strip limited 
by the dotted lines. 
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Fig. 23. The period-colour relation for stable unevolved sys- 
tems with Z = 0.02 (solid lines) and Z = 0.01 (dashed lines). 
The numbers denote the mass in solar units. 



stability problem. Recall that a rather large value for the ef- 
ficiency (/e = 0.1) is used. As shown in Table|3la reduced 
efficiency leads to an enhanced (but still small) degree of con- 
tact. 

3.2. Effects of metallicity changes 

The results obtained so far do not cover the range of observed 
systems. This can best be seen in the period-mass ratio dia- 
gram (Fig. [2}- Systems having large periods and intermediate 
or small mass ratios are absent. Stable systems with very small 
mass ratios are also absent. Since systems with these proper- 
ties are observed, effects neglected so far have to be taken into 
account. Here we consider the effects of metallicity changes in 
unevolved systems. 

The effects of a decrease in Z turned out to be small, apart 
from an increase in luminosity and effective temperature for 



3.66 
3.6 




19 P 

Fig. 24. The period-colour relation for stable evolved systems 
with Z = 0.02 and the masses M — IMq (dashed lines) and 
M = 1.8A/q (solid lines). The numbers denote the hydrogen 
content in the primary's centre. 

given mass and mass ratio. Figure^]shows examples in in the 
period-mass ratio diagram. Results for M = IMq represent 
very small masses where the difficulties are largest. Results for 
M = 1.8Mq represent typical masses. In both cases the influ- 
ence of Z is weak only. In particular, apparently the lower limit 
for the period depends only slightly on Z. The absence of sta- 
ble systems with large periods and intermediate or small mass 
ratios is even more pronounced when Z = 0.01. We conclude 
that metallicity effects cannot help to extend the range of stable 
systems in the period-mass ratio diagram. 

3.3. Effects of evolution 

Next we investigate evolutionary effects as the only remaining 
possibility to extend the range of stable systems. In the compo- 
nent i the simple hydrogen profile 



X 



0.7 



(0.7 - Xci) sin [kXi/ {2xc)] 



otherwise 



(12) 



will be assumed, depending on the central hydrogen content 
Xci and a parameter Xc- The standard value Xc = 0.5 will 
be used. Effects of changes in Xc will be investigated when 
analysing individual observed systems in Sect|3 

In the general case both components are evolved and the 
evolutionary status is described by two parameters Xd, Xc2 - In 
an approximate treatment evolutionary effects in the secondary 
can be neglected unless the mass ratio is large, and the evolu- 
tionary status is described by one parameter X^i- 

Concerning the choice of the circulation function we 
checked the invariance and found that the use of the standard 
circulation function is justified also in evolved systems. 

Evolutionary effects for systems with M ~ 1Mq,Z = 
0.02 are shown in Fig. |20l Except for the thin solid line all 
curves have been calculated neglecting evolutionary effects 
in the secondary. Heavy lines represent configurations with a 
given value for Xd- As before, solid/dashed lines represent 
configurations of positive/negative charge. Evolved configura- 
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tions with very large mass ratios {q ^ 0.8) have been omitted. 
With decreasing Xd the curves are shifting to the right, i.e. 
to larger periods, and the extent of stable configurations is in- 
creasing. Stable evolved systems with very small mass ratios 
exist for periods larger than about 0.3 days. 

On the dotted curves the angular momentum is constant. 
Since in the course of evolution Xd decreases and the angu- 
lar momentum does not increase, the period increases and the 
mass ratio decreases. The lower limit for the period found in 
unevolved systems is therefore valid also for evolved systems. 

In systems with large mass ratios the results are rough 
only since effects of nuclear evolution in the secondary are ne- 
glected. The influence of these effects can be seen if they are 
overestimated. This is the case in the thin solid line, which rep- 
resents stable configurations with Xci = Xc2 = 0.3 and q ^ 
0.8. The curves of constant angular momentum for Xd = Xc2 
almost coincide with the curves for Xc2 — 0.7. Accordingly, 
the main results obtained neglecting evolutionary effects in the 
secondary (increase of period and decrease of mass ratio in the 
course of evolution) remain valid if these effects are taken into 
account. 

A similar diagram for systems with AI = I.SAIq is shown 
in Fig.|^ Again it can be seen that the evolution (with or with- 
out loss of angular momentum) proceeds towards larger peri- 
ods and smaller mass ratios. Again stable evolved systems with 
very small mass ratios are found. Most important is the result 
that large periods and small mass ratios are possible among 
evolved systems. 

We conclude that the structure of contact binaries depends 
sensitively on evolutionary effects and that the range of ob- 
served systems can be explained only by these effects. 

3.4. The period-colour relation 

The period-colour relation of observed contact binaries (Eggen 
[196 1 1 967) provides an important observational test. In the 
period-colour diagram ("Figs. l22^24^ the observed systems are 
mainly found in the strip limited by the dotted lines. Figurel22l 
shows stable and unstable unevolved configurations with Z = 
0.02. Almost all configurations outside the strip are unstable. 

Let the strip be divided in two strips of the same width (up- 
per and lower strip, respectively). Figure |23l shows that stable 
unevolved systems cover roughly the upper strip, and Fig. |23 
shows that stable evolved systems cover almost the full strip. 
Accordingly, the full strip is covered by stable systems, and 
stable systems outside the strip are almost absent. 

These results lend support to the assumption of thermal 
equilibrium as well as to the treatment of the stability problem. 



4. Individual observational tests 

In this section individual observed systems with reliable data 
will be used as tests of the theory. First we shall ask whether 
such tests are severe and what we can learn from them. 
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Table 5. Models for BV Dra. 
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Table 6. Models for BW Dra. 



/e 






Z 




L 


Ri 


R2 


F 










5930 


1.40 


0.98 


0.55 


0.140 


10" 


-1 


200 


0.006 


5900 


1.33 


0.965 


0.537 


0.053 


10" 


'2 


200 


0.006 


5903 


1.35 


0.968 


0.541 


0.088 


10" 


-3 


200 


0.006 


5904 


1.37 


0.973 


0.546 


0.144 


10" 


-3 


200 


0.005 


5978 


1.44 


0.974 


0.547 


0.146 


10" 


-3 


100 


0.005 


5995 


1.46 


0.977 


0.550 
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4. 1. Observable quantities and degrees of freedom 

We begin with an observational difficulty. A light curve syn- 
thesis requires an assumption concerning the distribution of 
surface brightness. Usually a gravity brightening law is as- 
sumed (sometimes with different exponents in the two compo- 
nents) and the temperature difference is treated as an adjustable 
parameter (In particular, a negative temperature difference is 
adopted to reproduce the light curves of W-type systems.) The 
arbitrariness of this procedure is manifest. If a gravity brighten- 
ing law exists there is no freedom in the temperature difference. 

Actually there is no sound theoretical argument in favour of 
a correlation between local gravity and local effective tempera- 
ture. In a rotating single star without complications (e.g. spots) 
an (individual) gravity brightening law exists already for sym- 
metry reasons. In a contact binary symmetry arguments fail as 
well as physical arguments (c.f. Kahler & Fehlberg |l991> . and 
a detailed mapping of the surface appears to be necessary as 
recognized already by Hilditch et al. (■1988J . 
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Table 7. Models for 00 Aql. 
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In the present treatment the temperature difference is bound 
to be positive, and there is no freedom to explain the light 
curves of W-type systems with a negative difference. A nega- 
tive difference is also in conflict with infrared observations (e.g. 
Hrivnak ( 1989) and UV observations (e.g. Rucinski 19931. 

In view of these difficulties we consider the observed mean 
temperature as more reliable than the temperature difference, 
and the observed total luminosity L as more reliable than the 
luminosities Li, L2 of the components. 

Period P, mass ratio q, mass A/, mean temperature T^, lu- 
minosity L, radii Ri, and degree of contact F will be taken as 
observational quantities. They are not independent. Recall that 
Ri = XiA, where the separation A of the components is a func- 
tion of M and P on account of Kepler's law. The normalized 
radius Xi is a function of q and F. The luminosity L is known 
if radii and temperature are known. We are left with A^obs = 5 
independent observables P, q, M, T^, F. 

The degrees of freedom in a model are mass M, angular 
momentum (or period P), metallicity Z, central hydrogen con- 
tent in the primary Xd, and efficiency /e. (There is little free- 
dom in the temperature difference, and the effects of changes in 
a are unimportant). Accordingly, there are A^frcc = 5 degrees 
of freedom. 

Since A^froo > A^obs there is freedom to fit a model to 
the observations. The tests involving individual systems are 
therefore not able to provide evidence for thermal equilib- 
rium. Nevertheless they are important. Since A^frec = A'obs 
the model is well-determined. The structure of observed sys- 
tems can therefore be determined, assuming thermal equilib- 
rium, and the stability can be tested. In the case of instability 
either the treatment of the stability problem is insufficient or 
the system is in thermal imbalance. In the case of stability the 
efficiency can be calibrated. 

4.2. The system AB And 

For the system AB And (P = 0.331892 day) Hrivnak ( fT988l l 
determined the mass ratio q = 0.491, the mass AI — 1.5Mq, 
and otiier properties listed in the first line of Table |3 where 
temperatures are in Kelvin and luminosities and radii are in 
solar units. (This holds as well in the following tables.) The 
temperature is an average between the components. According 
to Hrivnak the uncertainty in L is large. Standard eiTors in the 
last digits are 1 for the radii and 3 for the degree of contact. 



0.6 - 




-0.6 - 

-0.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 

Fig. 25. The structure of AB And in the equatorial plane. 
Shown are the surface (solid line), the critical surface (dashed 
line), and the lower borders of the convective envelopes (dotted 
lines). The dots denote the centres of the components and the 
centre of mass. 

Treating the observed values for period, mass, and mass ra- 
tio as constraints and using Xd to adjust the mass ratio, a con- 
figuration is determined by the metallicity Z and the parame- 
ters /e, ATo. Models for AB And are summarized in Table |3 
The charge is positive and the models are indeed stable. This 
has been checked by evolutionary calculations following a 
small perturbations. The amplitude of the circulation is small 
(/ ~ 0.04). The models show that the theory is compatible 
with the observations. The results for the temperature suggests 
the metallicity Z = 0.011. The results for the degree of con- 
tact suggest a very small value for the efficiency (/e — 10~^). 
The model with these two values and with ATc = 200 K is in 
excellent agreement with the observations, except for the radii 
which are somewhat too small. We shall return to this point. 

The geometry of the system in the equatorial plane is shown 
in Fig.|25] Note that the fractional extent in radius of the con- 
vective envelope is much larger in the primary than in the sec- 
ondary. The fractional extent in mass is also much larger. 

We checked also the effects of changes in the chemical pro- 
file. In the models listed in Table |3 calculated with the stan- 
dard profile (xc = 0.5), the hydrogen content Xd is somewhat 
lower than 0.4. Models with Xc ~ 0.3 give similar results with 
Xd — 0.2. Accordingly, the choice of the chemical profile is of 
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Fig. 26. log P as a fuction of the potential in the common en- 
velope of system 2. The strong solid/dashed line represents the 
primary/secondary. The dotted line represents a modified run 
of pressure in the secondary (see text). 

minor importance. The amount of hydrogen already converted 
to helium is also similar and can be used to estimate the age of 
the system. 

4.3. The system BV Dra 

The contact binaries BV Dra and BW Dra are particularly in- 
teresting since they form a visual binary. Precise observational 
results have been obtained by Kaluzny & Rucinski (119861 . We 
begin with BV Dra. Period (0.350067 day), mass (1.47M0) 
and mass ratio {q = 0.411) are again treated as constraints. 
Other observational results are listed in the first line of Table|5] 
Standard errors in the last digits are 21 for the luminosity, 1 for 
the radii, and 27 for the degree of contact. 

The models of BV turned again out to be stable. Comparing 
the temperatures with the observed value we obtain the metal- 
licity Z = 0.006, in accordance with the result of Kaluzny & 
Rucinski. The degree of contact suggests again a very low ef- 
ficiency (/e — 10^"^). The model with these values and with 
ATc = 200 K is in very good agreement with the observations, 
again apart from the fact that the radii are somewhat too small. 

4.4. The system BW Dra 

The constraints for this system are P = 0.292167 day, M — 
1.18Mq, q = 0.280. Other observational results are listed in 
the first line of Table |6l Standard errors in the last digits are 
1 for the radii, 15 for the luminosity, and 34 for the degree 
of contact. Since no mean temperature has been reported by 
Kaluzny & Rucinski, a temperature somewhat smaller than the 
primary's polar temperature (5980 K) has been adopted, as it is 
the case in BV Dra. 

The models are again stable. This has been checked in evo- 
lutionary calculations. The temperature suggests a metallicity 
Z = 0.005 . . . 0.006 which is very similar to the metallicity of 
BV Dra. This is as expected. The degree of contact suggests 
again the efficiency /e — 10^"^. Adopting these values and a 



temperature difference ATe = 200 K, the agreement between 
theory and observations is very close. 



4.5. The system 00 Aql 

The system OO Aql, observed by Hrivnak ( I1989t . has a very 
large mass ratio (q — 0.843). Other constraints are M = 
1.92Mq, P = 0.506789 day. The observed temperature dif- 
ference ATc = 65 K has been added as a constraint. Other 
observed properties are listed in the first line of Table0 

Since q is large, evolutionary effects in the secondary are 
important. (If they are ignored, the constraints require a neg- 
ative amplitude /.) Table contains properties of configura- 
tions for different combinations of Xd and Xc2- These com- 
binations are approximately determined by the condition that 
/ is positive and small. Note that Xc2 is only slightly larger 
than Xci and that these values depend on the choice of Xc- 
Evolutionary effects are large in this system, in accordance 
with the result of Hrivnak. The configurations in the last two 
lines are in excellent agreement with the observations. The 
radii are slightly smaller than those given by Hrivnak. The 
agreement is nevertheless perfect. Hrivnak's radii are volume 
radii (c.f. Mochnacki 1984 1 and thus larger (in the present 
case by about 1 .4 percent) than the radii in a spherically aver- 
aged treatment of the components (c.f. Kahler 1986 1. (A sim- 
ilar difference in radii occurs in the systems discussed above). 
The metallicity (Z = 0.013. ..0.014) turns out to be larger 
than estimated by Hrivnak. The efficiency is extremely small 
(/e ~ 10^^), and the amplitude is rather large (/ ~ 0.09). 

Again the stability was tested in evolutionary calculations. 
The two best models (in the last two lines of Table^} are stable, 
but the other models are unstable. Apparently an extremely low 
efficiency (or a rather large degree of contact) is necessary for 
stability. This is suggested also by other models not listed in 
TableQ 



4.6. Conclusions 

Several observed systems with well-determined parameters 
have been used as tests of the theory. In all cases an excellent 
agreement between theory and observations is obtained. This is 
as expected from the number of free parameters. In all systems 
evolutionary effects are important. 

All systems turned out to be stable. This lends support not 
only to the transport equation used in the present model but also 
to the assumption of thermal equilibrium and to the treatment 
of the stability problem. 

The efficiency was determined to be very small (/e — 
lO^'^ . . . 10^^). This shows that the velocities of the internal 
mass motions in the neighbourhood of the inner Lagrangian 
point are much smaller than the sound velocity. 
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Fig. 27. Surfaces of constant pressure (solid line), constant tem- 
perature (dashed line), constant entropy (dash-dotted line), and 
constant density (dotted line) in a radiative region of a rotating 
star. 

5. Towards an understanding of the internal 
mass motions 

5. 1. Internal mass motions In the common envelope 

So far Roche geometry has been assumed and the effects of in- 
ternal mass motions on the equation of hydrostatic equilibrium 
have been neglected. This is justified as a close approximation 
to real systems, but a closer look into the layers in the com- 
mon envelope reveals an inconsistency. This inconsisteny can 
be removed by a slight modification, and this procedure offers 
insight into the internal mass motions in the common envelope. 

As noted in Sect. 12.51 in the present approximate treatment 
the pressure difference AP on equipotential surfaces is nega- 
tive apart from the surface. For the layers in the common en- 
velope of system 2 this is shown by the strong lines in Fig. 1261 
Since the large-scale mass motions are driven by pressure dif- 
ferences on equipotential surfaces, and since this difference is 
negative in each layer, mass is flowing from the secondary to 
the primary but not from the primary to the secondary. This is 
incompatible with a configuration in thermal equilibrium with 
a vanishing net rate of mass transfer. 

Can this inconsistency be removed by a slight modification 
of the structure equations, taking into account internal mass 
motions? The pressure difference in the common envelope is 
largest (apart from the sign) in the layers just above the critical 
surface, causing a mass flow from the secondary to the primary 
in these layers. In a slightly changed and more consistent treat- 
ment this feature is preserved. We need therefore a flow from 
the primary to the secondary in the outermost layers in order to 
obtain a vanishing net rate. The resulting velocity field is of the 
type discussed by Nariai ( 1976j and Zhou & Leung ( 1990 ). 

The effects of this velocity field on the run of pressure in 
the common envelope are determined by a modified surface 
condition (the generalization of the Roche equipotential condi- 
tion), a change of the effective gravity, and (since the motions 
are expected to be turbulent) the presence of a turbulent pres- 
sure. We shall discuss these effects in turn. Note that the Roche 



potential is well defined and that the difference between Roche 
potential and effective potential is negligible in the common 
envelope. The pressure can therefore be treated as a function 
of the potential. Since we ask for slight changes, we may as- 
sume that the change of a curve in Fig.|26lis a superposition 
of a horizontal shift (caused by the change in the surface con- 
dition) and a change in the pressure gradient (caused by the 
change in the effective gravity and by the turbulent pressure). 
To reduce the inconsistency we need a shift to the left/right or a 
decrease/increase in the gradient of the secondary's/primary's 
curve. 

Since the velocity field has a reversing layer, from 
Bernoulli's equation the surface condition 



Ml 

2 



^2 + 



Ml 

2 



(13) 



is obtained, where Vi is the velocity at the surface of the com- 
ponent i and the angle brackets denote a spherical average, cf. 
Kahler (11995V Since the pressure vanishes at the surfaces, the 
surface condition acts as a boundary condition for the run of 
pressure. Assuming a horizontal shift of the curves, it can be 
veryfied that a reduction of the inconsistency requires that 



(14) 



i.e. that the velocities in the secondary are larger than in the pri- 
mary. The dotted line in Fig.|26lshows an example for the result 
of a horizontal shift of the secondary's curve. It is manifest that 
the inconcistency can be removed in this way. 

Consider next the effective gravity. In the presence of inter- 
nal mass motions gi has to be replaced by 



Gr 



9i 



with il? > uj"^. 



(15) 



i.e. the effective gravity is reduced (Kahler fl995> . This implies 
that also the pressure gradient is reduced. If this occurs in the 
secondary, again the inconsistency is reduced since the pressure 
at the critical surface decreases. 

It remains to discuss the turbulent pressure. Let Pg i be the 
gas pressure and ~ fiPg.i the total pressure. Neglecting 
the radiation pressure we have = 1 in the absence of in- 
ternal mass motions and ~ 2 if the turbulent velocities are 
comparable with the sound velocity. Combining the equation 
of hydrostatic equilibrium and the equation of state we find 



dr,; 



InP 



(16) 



Accordingly, an increase in /; has the same effect as a decrease 
in the effective gravity. Turbulent pressure reduces the pressure 
gradient. If this occurs mainly in the secondary, the inconsis- 
tency is again reduced. 

The discussion of the three effects has shown that the in- 
consistency can be removed by mass motions which are larger 
in the secondary than in the primary. Pressure differences in the 
layers avove the critical surface provide the driving mechanism 
for the velocity field in the common envelope. The result is a 
field with a reversing layer, with motions from the secondary 
to the primary in the region above the critical surface and from 
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Fig. 28. The normalized gravity 7 at the top of the radiative 
interior in the secondary of AB And, in dependence of the an- 
gular distance (in degrees) i? from the x-axis. Shown is the run 
of 7 in the plane containing the rotation axis (solid line) and 
in the equatorial plane (dashed line). The dotted line shows the 
averaged gravity. 



the primary to the secondary in the surface layers, and with ve- 
locities which are larger in the secondary than in the primary. 
Martin & Davey (|1995) obtained in hydrodynamic simulations 
a velocity field with the same properties. 

The arguments presented in this subsection concern late- 
type as well as early-type systems. They are not conclusive 
since in a spherically averaged treatment of the components 
degrees of freedom of the internal mass motions escape notice. 

5.2. Circulation currents in single stars 

Consider first the radiative envelope of an early-type single star 
in almost uniform rotation. Near the boundary to the convective 
core as well as near the surface viscous effects are important 
as discussed by Tassoul & Tassoul (1995 1, but in the bulk of 
the envelope the original inviscid solution of Sweet ( 1950 1 is a 
reasonable approximation. Accordingly, on an isobaric surface 
the temperature is higher around the poles (where the gravity is 
larger) than around the equator (where the gravity is smaller). 
The resulting topology of the surfaces of constant thermody- 
namic quantities is shown in Fig.|^ The eccentricity is small- 
est in the surface of constant temperature and largest in the sur- 
face of constant entropy. We shall need the result that also on 
a surface of constant entropy the temperature is larger around 
the poles than around the equator 

Let l{s) be the luminosity carried through the surface s ~ 
const, which is the sum 



lis) = /rad(s) + /circ(s) 



(17) 



of the radiative luminosity ^rad and the circulation luminosity 
Zcirc- In thermal equilibrium and outside the core with nuclear 
burning we have dl/ds — and thus 



Assuming stationary currents the heat equation reduces to 

divF = -pTu-grads, (19) 

where F denotes the radiative flux and u is the circulation ve- 
locity. We integrate this equation over the volume enclosed by 
the surfaces s = si and s — S2, where the difference S2 — si is 
positive and infinitesimal. The result is 



^rad(s2) - 'rad(si) = -(s2 - Si) 



pT (u • n) dA, (20) 



where n = grads/|grads| is the outward directed normal on 
the surface s — S2 and dA is an element of this surface. Let T 
be the average of T on this surface and AT = T — T. Since 
the net mass flow trough the surface vanishes we obtain 



dl 



rad 



j ATp (u • n) dA. 



(21) 



Since AT and u • n have the same sign almost everywhere, 
positive around the pole and negative around the equator, the 
integrand is positive almost everywhere. Therefore d^rad / ds is 
negative, and from Eq. ( I18t we see that dZcirc/ds is positive. A 
slightly generalized discussion shows that this result remains 
valid in the presence of nuclear burning. 

We conclude that radiative layers are sources of the circu- 
lation luminosity. The circulation luminosity is positive and a 
monotonically increasing function of s or of the mass variable. 
These results concern the bulk of the radiative envelope where 
viscosity is unimportant, but not the regions at the borders of 
the envelope. Recall that sinks near the surface are needed. 

5.3. Circulation currents in contact binaries 

First let us ask for the changes of the gravity g on a level sur- 
face in the secondary's interior The mass TO2 enclosed by this 
surface is known from the spherically averaged treatment of the 
system. In a close approximation g is (apart from the sign) the 
gradient of the potential 

GMi Gma 1 



dl 



(22) 



d/circ/ds = -d^iad/ds. 



(18) 



where di is the distance from the centre of the component i and 
d is the distance from the rotation axis. 

Taking the secondary of AB And as an example. Fig. |28l 
shows the normaUzed gravity 7 = gA^ / (GM) at the bottom of 
the convective envelope, i.e. at the top of the radiative interior. 
The dashed line shows 7 in the equatorial plane in dependence 
on the angular distance •& from the .T-axis. (As usual, the x-axis 
is taken to lie along the line of centres of the components, with 
the origin in the primary's centre.) The solid line shows the run 
of 7 in the plane containing the rotation axis, and the dotted line 
denotes the averaged gravity obtained in a spherical treatment 
of the components. The curves show that the gravity changes 
almost by a factor of two. On level surfaces deep in the interior 
the variation is smaller, but on the surface having (in a spherical 
treatment) the fractional radius 0.5 the gravity still changes by 
about 10 percent. 

On each level surface the gravity is largest at the poles 
(solid curve for § = 90°), small near § = and smallest 
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when I? = 180°, i.e. when closest to the inner Lagrangian point. 
On account of tidal synchronization the rotation is known to be 
almost uniform. Although the departures from sphericity in the 
outer layers are large and although the geometry is more com- 
plex than in single stars, the essential features are similar. We 
therefore expect rising motions at the poles and sinking mo- 
tions at the equator, at least near the x-axis. Similar motions 
are expected also in the primary. 

This is supported by results on radiative regions in uniform 
rotation (Mestel .l966j . They predict a double-cell pattern with 
a circulation inversion on the level surface with the density 
p = uj^ I (27rG). In AB And the outer cell is absent in the pri- 
mary (since the critical density occurs in the convective enve- 
lope) and unimportant in the secondary since the bulk of the ra- 
diative interior is covered by the inner cell. In this cell the circu- 
lation is upward/downward in regions of higher/lower gravity. 
Mestel's result concerns an idealized problem since barotropy 
is assumed and viscosity is neglected. Questions about the va- 
lidity of Mestel's result (Tassoul & Tassoul .l995il concern the 
outer cell and the region around the circulation inversion, but 
not the bulk of the radiative region where Sweet's (I1950> re- 
sults are confirmed. 

Accordingly, it is certain that on on an isobaric surface the 
density is lower in regions of outward motions than in regions 
of inward motions. Therefore the discussion of the preceding 
subsection applies, showing that there are sources of the circu- 
lation luminosity throughout the bulk of the radiative interior. 
This is in line with the basic assumption made in Sect[2 It is 
also in line with the result from Sect. 12.21 that stable systems 
have sources not only in the outer layers but also deep in the 
interior. The corresponding sinks occur either in the outermost 
radiative layers or in the common turbulent envelope. 

6. Summarizing and concluding remarl<s 

The present investigation was based on the assumption that the 
circulation luminosity in the secondary's interior is positive and 
large enough to enable thermal equilibrium. The first part of 
this assumption can be abandoned. As shown in the preced- 
ing section the circulation luminosity is certainly positive in 
the radiative interior, not only in the secondary but also in the 
primary. In systems with intermediate or small mass ratios the 
effects in the primary are small while the effects in the sec- 
ondary are important since the core luminosity is small. In sys- 
tems with large mass ratios (e.g. OO Aql) the neglect of the 
circulation luminosity in the primary introduces a slight error. 

The resulting treatment of contact binaries is simple. 
Sources of the circulation luminosity deep in the secondary's 
interior are necessary, but details are unimportant. Free param- 
eters are the temperature difference between the components 
ATc, bound to be positive and small, and the efficiency /e of 
the energy transfer between the components. The fractional ex- 
tent of radiative regions is larger in the secondary than in the 
primary. In the course of evolution the period increases and the 
mass ratio decreases. 

Adopting standard values for the parameters, a survey of 
unevolved and evolved contact configuration has been ob- 
tained. The results are successful in explaining the basic ob- 



servational facts, in particular the preference for shallow con- 
tact, the minimum period and the minimum temperature, the 
existence of stable configurations with mass ratios very close 
to unity, and the period-colour relation. Since stability consid- 
erations are essential in these observational tests we consider 
the results as strong arguments not only in favour of thermal 
equilibrium but also in favour of the treatment of the stability 
problem. 

Assuming thermal equilibrium, models of individual ob- 
served systems with reliable data are well-determined apart 
from some freedom in ATo. When stable they can be used 
to to calibrate the efficiency and to determine metallicity and 
age. The models obtained so far are stable, which lends again 
support to the assumption of thermal equilibrium. Evolutionary 
effects are important. The efficiency is very small (/e — 

Concerning the hydrodynamic problem, arguments have 
been presented that the velocity field in the common envelope 
has a reversing layer with motions from the secondary to the 
primary in the layers just above the critical surface and from 
the primary to the secondary in the surface layers. 

Details of the circulation, in particular the circulation func- 
tion, can be determined only in a three-dimensional discussion. 
In comparison with rotating stars the problem is more diffi- 
cult since the geometry is more complex. In other respects the 
problem is simpler since rotation law and geometry are known 
in a close approximation. This problem will be investigated in 
a forthcoming paper, again assuming thermal equilibrium. 

Although we have strong arguments in favour of thermal 
equilibrium we have no proof, and it is uncertain whether a 
proof is possible. 
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